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Abstract 

A method of reconstruction of the top quarks produced in the process e + e~ — > 
it — ► 6 jets at a Linear Collider (LC) is proposed. The approach does not involve 
a kinematic fit, as well as assumptions on the invariant masses of the dijets 
originating from the decays of W bosons and, therefore, the method is expected 
to be less sensitive to theoretical and experimental uncertainties on the top- 
mass measurement than traditional reconstruction methods. For the first time, 
the reconstruction of the top quarks was investigated using the full LC detector 
simulation after taking into account the background arising from QCD multi-jet 
production. 



1 Introduction 

The standard procedure of the top-quark measurements in the process e + e~ — > it — > 
^<?i<?2<?3<?4 — > 6 jets takes advantage of a fully reconstructed final state. The decay sig- 
nature is characterized by the production of six hadronic jets, therefore, a background 
arising from the standard six-jet QCD processes is expected. Nevertheless, since the 
all-hadronic top decay has the largest branching ratio (~ 44% of all it decays), it is con- 
sidered as one of the most promising at the TEVATRON, LHC and a Linear Collider 
(LC). 



A reduction of the background from QCD multi-jet events can be achieved by 
identification of ^-initialized jets (the so-called b— tagging). Furthermore, one can use 
the requirement that the invariant mass of two jets (not associated with the b quarks) 
is consistent with the decay of the W bosons. These two methods were used for the 
top reconstructions at the TEVATRON (see, for example, pQ), and also are considered 
for the LHC and LC experiments 0HIEI- 

The reconstruction of the top quarks in the all-hadronic decays is known to be 
affected by many theoretical uncertainties: the extraction of the top-quark pole mass 
has a theoretical uncertainty of around 300 MeV and cannot be determined with a 
precision better than O(Aqcd) [§]• I* 1 addition, incomplete knowledge of the hadronic 
final state leads to an uncertainty of a few hundreds MeV (3). For the latter, one 
of the most significant uncertainties is related to the reconstruction of the W mass, 
My/, after minimizing | Mjj — M\y |, where Mjj is the dijet invariant mass. The way 
how the large Breit-Wigner tails of the decay W — > 2 jets is treated when the top 
quarks are reconstructed directly affects the reconstructed top mass [5]. A significant 
impact on the reconstructed top mass may come from non-perturbative phenomena, 
such as colour rearrangement j7j and Bose-Einstein effects jH], which are expected to 
shift the reconstructed masses in the W + W~ and ti production. Thus, any approach 
involving the reconstruction of W — > 2 jets is bound to lead to a systematic shift for the 
top reconstruction. For example, it has been shown that the present implementation 
of the Bose-Einstein effect in a Monte Carlo model can contribute to a shift in the 
reconstructed top mass through a distortion of the Mjj spectrum 0. Similar effects are 
expected for the color-reconnection phenomenon which is anticipated to be significant 
for the LC energies. 

In this article, a simple reconstruction method which does not involve the direct 
measurement of the decay W — > 2 jets is proposed. This method is a very general, 
and is suitable for the reconstruction of any process in which a particle (V) and an 
anti-particle V are produced and then decay as: 

e + e - _> V + v — ► v a V a + v b V b 

V a — >Vi + v 2 , v b — >v 3 + v 4 , 

without a prior knowledge on the masses of the intermediate particles, denoted as V a 
and Vb. The effectiveness of this method is based on the assumptions that the initial- 
state particles, V and V, have similar (but not equal!) masses, | My — My |<< M y( y), 
and that the simplicity of e + e~ annihilation allows to use the momentum conservation, 
Py + P v = 0. These two requirements, together with the fact that both initial particles 
decay similarly into three other particles, are essential in the reconstruction of the 
invariant masses of V and V through six jets in the final state. 

2 Reconstruction procedure 

We consider e + e~ collisions in the laboratory frame. As a first step, six jets in the event 
have to be reconstructed. We use the k± (Durham) algorithm [D], requiring exactly six 
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jets in every event, without any specific cuts on the transverse momenta and rapidity 
of the jets. 



2.1 Event selection 

To preselect hadronic events, we use the typical LEP cuts which take advantage of the 
energy-momentum conservation requirements: 

\E vis I Ejjjjj I a I Y,PTi I < a m 

I —=■ — L \< 1\ E , _» < Ap L , _> < IXpT, {i-) 

V s 1^ I Pi I 2-J -Pi I 

where P„j S is the visible energy, (pYi) is the longitudinal (transverse) component of 
momentum of a final-state particle. The Ap, Ap^ and App are small adjustable param- 
eters. The top events are characterized by a large amount of missing energy /momentum 
due to undetected neutrinos, therefore, the selection (JTJ) is especially important as it 
helps to reject events with a significant fraction of neutrinos [5]. 

The jets in the all-hadronic top-decay channel should be well separated in the 
transverse momenta. The best way to do this is to apply a restriction on the values 
of ?/g ut used in the Durham algorithm to resolve six jets. We accept only events if 
yl ut > A y , with A y being a parameter to be found using a Monte Carlo simulation. A 
similar cut was applied for the JADE jet algorithm in the study of the all-hadronic top 
decay using the LCD Fast Simulation [3]. 



2.2 Top reconstruction 

We start with the initial list of six jets with momenta pi. These jets are merged 
into groups, with three jets in each group. Thus, every event can be considered to 
consist of pairs of the three-jet groups with momenta {P/(l), P/(2)}, where Pr(l) = 
Pi,j,m — Pi + Pj + Pm 1S the four-momentum of each group, and P/(2) was obtained 
analogously using the rest of the jets in the same event. In total, there are twenty 
three-jet groups, which can be arranged into ten three-jet pairs with non-identical 
momenta, i.e. / = 1, ... 10. Ideally, from these ten pairs of three-jet clusters, one 
should accept only one three-jet pair corresponding to two initially produced particles. 
In practice, it is difficult to find the correct assignments of jets due to fragmentation 
and detector effects. Therefore, a special selection should be performed in order to 
minimize the number of three-jet pairs. Assuming that the produced particles have 
similar masses, we accept such three-jet pairs if the invariant masses of the three-jet 
groups satisfy the condition |M/(1) — Mj(2)| < Am, where Mj is the invariant mass 
of a three-jet group with the momentum Pj, and Am is a free parameter constraining 
the invariant masses of three jets inside each pair. 

Next, energy-momentum conservation can be used to set an additional constraint 
on the momenta of three-jet groups within each three-jet pair. We accept only three-jet 
groups which are produced back-to-back, i.e. we require | P/(l) +P/(2) |< Ap, where 
Ap is again an adjustable parameter. 

The following additional selection criteria may also be used: 
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• The list of three-jet groups contains top-quark candidates, which decay as t — > 
bW — > bqq. In order to reduce the number of possible three-jet pairs, one can 
use an efficient double fc-tagging, requiring one and only one 6-initiated jet within 
each three-jet group. All three-jet pairs which have at least one three-jet group 
with two 6-initiated jets have to be removed. This leads to six possible three-jet 
pairs, compared to ten if no 6-tagging is used. 

• Each three-jet group should contain two jets coming from the decay of the W 
boson. Thus, one can remove such three-jet groups which cannot be associated 
with the VT-decay hypothesis. This can be done by accepting such three-jet pairs 
which have the invariant masses of any two jets inside the region | Mjj — My/ |< 
A y/. This requirement is not advocated in this article as it makes the specific 
physics assumptions on the invariant mass and width of the W boson. 

These two additional criteria are optional; in fact, we will show that even without 
them one can obtain a robust reconstruction procedure. 

2.3 Monte Carlo studies 

To illustrate the method outlined above, we first apply it to the all-hadronic top decays 
in e + e~ annihilations at the centre-of-mass energy of y/s = 500 GeV generated using 
the PYTHIA 6.2 model [TP] . The contribution from the initial-state radiation (ISR) 
can be rather significant at LC energies, thus this effect was included in the simulation. 
The mass and the Breit-Wigner width of the top quarks were set to the defaults values, 
175 GeV and 1.39 GeV, respectively. The particles with the lifetime more than 3 cm 
are considered to be stable. Neutrinos were removed from the consideration. 

In order to reconstruct the top quarks using the proposed method, one should define 
the following: a) the parameters for the event selection: Ap, Apl, Apr, A y ; b) the 
parameters for the final reconstruction: Am and Ap. The event-selection parameters 
should not be large in order to insure small contaminations from neutrinos resulting 
to a broad and asymmetric Breit-Wigner peak for the three-jet mass distribution. 
In this study, we will apply rather tight cuts: Ap = Ap^ = App = 0.02. The 
parameter A y requires an additional study: Fig. ^ shows the distribution for the all- 
hadronic top decays and for e + e~ inclusive events (but without the e + e~ — > tt — > 6 jets 
process) generated with the PYTHIA model. The top events are characterized by 
ycut > = 2 • 10~ 4 , thus this value is used for the studies below. 

The parameter Am controls the maximum difference between masses of two top- 
quark candidates. It was set to 30 GeV, which is large enough to insure a reasonable 
event acceptance, but sufficiently small to exclude wrong assignments of jets. The 
parameter Ap = 5 GeV further reduces the contributions from unmeasured neutrinos. 
Note that the latter cut is not very important for the generator-level study, since the 
event selection is already very strong. In practice, such cuts should not be very tight 
to obtain a reasonable acceptance after a detector reconstruction. 

Figure 121 shows the three-jet invariant masses after the selection procedure described 
above, but without the requirements of double 6-tagging, and also without the cut on 
| Mjj — Mw | • Two different functions were used to fit the peak: a) the Breit-Wigner 
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distribution (with the fixed width of 1.39 GeV) convoluted with a Gaussian and b) the 
Breit-Wigner function alone. The first-order polynomial is used in order to describe the 
background in both cases. The motivation for the Gaussian convoluted with the Breit- 
Wigner function comes from the assumption that contributions of many independent 
effects (parton-shower splittings, hadronization and resonance decays) should yield a 
Gaussian-like distribution. As seen from Fig. |2l the best fit can be obtained by using 
the Breit-Wigner function alone, while the Gaussian does not describe adequately the 
tails of the mass distribution. Based on this empirical observation, we will use the 
Breit-Wigner function to fit the mass spectra before the LC detector simulation. 

The fact that the Gaussian distribution fails to describe the mass spectrum might 
indicate a strong correlation between different effects leading to non-Gaussian broad- 
ening of the natural width. It was verified that a sum of two Gaussian distributions 
gives a better fit than the Gaussian distribution alone, but the Breit-Wigner still gives 
the best fit. It is likely that a sum of many Gaussian distributions should be used to 
fit the mass spectrum. In this case, an unknown weight for each Gaussian contribution 
would reflect the fraction of events with a particular property. 

One interesting feature of the invariant-mass distribution shown in Fig. |2] is that it 
has no background in the region of small three-jet masses, where the major contribu- 
tion from the QCD multi-jet background is expected. This can be explained as follow- 
ing: any misassignment of jets at a centre-of-mass energy far from the top-production 
threshold usually increases the observed three-jet invariant mass, since in this case the 
misassignment jet usually has a significant angle with respect to the rest of the jets in 
the three-jet group. This ultimately increases the three-jet invariant masses. 

As was mentioned in the introduction, this method is expected to be less sensitive 
to the hadronic-final state uncertainties which usually arise when the invariant mass 
of the dijets from the decay W — > 2 jets is reconstructed. Fig. |U shows the three-jet 
invariant-mass distributions without and with the Bose-Einstein interference included. 
We use the so-called "BE32" type of the Bose-Einstein simulation, which is included in 
the PYTHIA 6.2 model ^U]- It is clear that no any significant shift for the reconstructed 
top mass can be attributed to the Bose-Einstein effect when the proposed method is 
used for the top reconstruction. 

3 Top reconstruction using the full 
detector simulation 

The reconstruction of the top quarks after the full detector simulation was performed 
in a few steps: 

1) Fully inclusive e + e~ annihilation events were generated with the PYTHIA 6.2 
model, which was based on the default parameters. Events with the W + W~ Z and ZZZ 
production leading to six jets in the final state were not included in the simulation, as 
they have a small contribution (~ 2% — 4%). In total, 270k events were generated. 
This sample approximately corresponds to twenty days of the LC data taking. 

2) From this sample, e + e~ — > ti — > everything events were passed through the 
full GEANT-based detector simulations. The number of such events was 9.6k. The 
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TESLA detector simulation based on the BRAHMS 305 program was used 
for the event reconstruction. Of particular importance for the present study are the 
tracking system based on the Time Projection Chamber (TPC) and the calorimeter 
(CAL). The TPC is surrounded by the CAL, which is longitudinally segmented into 
electromagnetic and hadronic sections. The electromagnetic calorimeter is based on 
tungsten absorber and silicon diode pads, while the hadronic part is an Fe/scintillating 
tile calorimeter. All particles in events were reconstructed using a combination of track 
and calorimeter information that optimizes the resolution of the reconstructed jets [T3] . 

3) The rest of the sample (without the top quarks) was treated differently: Energies 
of the final-state particles in such events were smeared around the true values accord- 
ing to the Gaussian resolution functions in order to imitate the detector response. 
Charged hadrons (with energy denoted by E c h) were smeared using the Gaussian dis- 
tribution with a = 10~ A E^. h . Analogously, photon energies (E^) were smeared using 
a = 0.12-i/E^, while the energies (E n ) of neutral hadron were distributed according 
to a = OAy/E^. Particles close to the beam pipe were removed, to reproduce the 
geometrical acceptance of the LC detector. 

The selection and the reconstruction of top quarks at y/s = 500 GeV are based on 
the parameters given in Tabled). These parameters are the same for the generated and 
reconstructed event samples. For the former sample, neutrinos were removed, as well as 
particles which are inside the beam pipe. The cuts on the energy/momentum imbalance 
were determined after the study of the detector resolution, and were optimized to obtain 
a reasonable efficiency of the reconstruction. Figure El illustrates the distributions of 
the total event energy before and after the detector simulation. A significant missing 
energy due to neutrinos and the ISR is observed even before the event reconstruction. 
The event-selection cuts are by a factor two-three larger than the resolution on the 
corresponding variable. Events were removed if at least one lepton is found with an 
energy above 20 GeV; the motivation for this cut is clear from Fig. |3J 

Figure shows the three-jet invariant masses for the reconstructed all-hadronic 
decays in the process e + e~ — > it before the detector simulation. The width of the 
Breit-Wigner distribution, 9.5 GeV, is somewhat larger than that shown in Fig- El since 
now a larger contribution from the missing energy/momentum is allowed to obtain a 
higher selection efficiency This width will be used in the following figures for the 
Breit-Wigner function. 

Figure El shows the three-jet invariant masses after the full detector simulation and 
event reconstruction. The background events from the non-top continuum were added 
after the Gaussian smearing. The cuts are the same as for Fig. |3J except for the cut 
on ?/g ni . A significant reduction of the background can be obtained by selecting events 
with jjQ Ut > 2 ■ 10~ 4 , as seen from Fig. Note that this cut is the same as for the 
generator-level studies, since the LC detector simulation does not distort significantly 
the yg ut distribution. 

Additional assumptions can further increase the signal-over-background ratio. For 
example, Fig. |H] shows that the cut on the invariant mass of two jets inside each three- 
jet group, | Mjj — Mw |< 15 GeV, reduces the background. Further, a check using 
the generator-level reconstruction indicated that an efficient double ^-tagging can de- 
crease the background by ~ 50%. This, however, has no a significant impact on the 
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Gaussian width, which is 5.0 — 5.5 GeV for the the current detector design and event 
reconstruction. 

4 Top reconstruction at = 800 GeV 

The proposed method is well suited for a higher centre-of-mass energy of e + e~ col- 
lisions. The reconstruction is expected to be even more reliable: Since the hadronic 
jets are better collimated along the momenta of the produced top quarks, any wrong 
jet assignment significantly increases the invariant mass for such a miss-reconstructed 
top candidate. This should lead to a large invariant- mass difference between miss- 
reconstructed three-jet groups, therefore, such groups which will be removed by the 
cut on |Mf(l) — M/(2)| more effectively. 

At higher energies, however, missing event energy /momentum is larger due to higher 
energies of neutrinos. Therefore, the selection cuts should further be optimized to 
obtain a reasonable selection efficiency. Table ^>) shows the event-selection and top- 
reconstruction cuts applied for \/s = 800 GeV. We also use the same reconstruction 
cuts as for yfs = 500 GeV to be able to compare the results with the previous founding. 
Note that the y^ 1 distribution for = 800 GeV is somewhat different than that shown 
in Fig. ^ now the distribution extends down to ?/g n * = 10 -5 , and has a small peak at 
yc ut ~ 10~ 4 . This is due to a better collimation of jets originating from the top quark. 
It was found, however, that a decrease of A y for events at yfs = 800 GeV has a small 
impact on the final selection efficiency. 

Figure shows the three-jet invariant masses for e + e~ — > ti events generated with 
PYTHIA before the detector simulation. The PYTHIA event sample has the same 
size as for Sect. El The signal-over-background ratio is significantly larger for events at 
a/s = 800 GeV than for ^Js = 500, as expected. 

Note that the efficiency of the top selection is less by a factor two than for e + e~ 
annihilation at -y/i = 500 GeV. It was found that most of the events were rejected by 
the cut |-P/(1) + P/(2) | < Ap, thus this cut is the main in reducing the fraction of 
events with significant missing energies from neutrinos. A check without the cut on the 
energy/momentum imbalance shows a very similar result as for the case without these 
cuts included. Generally speaking, the cuts on the event imbalance can be applied 
mainly to reduce the computational time before the top reconstruction, rather than to 
improve the signal-over-background ratio. 

Figure El shows the reconstructed signal after the detector simulation and when 
events from the continuum were included. No any cuts on the dijet invariant mass 
were applied, i.e. this figure is equivalent to Fig. [7| for the lower centre-of-mass energy. 
The signal for = 800 GeV has significantly less background, however, the number 
of the reconstructed top candidates is also smaller due rejection of events with large 
energy losses. 

It is important to notice that the width of the Gaussian distribution is larger for 
y/s = 800 GeV than for -^/i = 500 GeV. There are a few reasons for this. First, the 
energy resolution of the CAL is worse for higher jet energies. This, however, cannot 
completely explain the observed increase in the mass width. Again, let us remind that 
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jets from the decay of a single top are better collimated at s/s = 800 GeV, thus a larger 
degree of overlaps between energy deposits in the calorimeter is expected that makes 
more difficult to reconstruct the energy-flow objects. Furthermore, one should expect 
a stronger leakage for neural particles outside the calorimeter. 

5 Snagging the top quarks using a neural network 

Neural networks (NN) have seen an explosion of interest over the last years, and have 
been successfully applied for paten recognition problems in particle physics. Based on 
the top-reconstruction approach described above, in this section we will devise a NN 
method to select top events using the hadronic-final-state signatures of the all-hadronic 
decays. 

As a first step, one should define the variables for NN input which are likely to 
be influential. The event-shape variables represent the most efficient way to separate 
the ti from the continuum [THE]. The thrust (T) [T^j, major (Mj) and oblateness 
(Ob) event characteristics will be our primary choice. The allowed range for thrust 
is 0.5<T<1, such that the isotropic events have T ~ 0.5, whilst a 2-jet configuration 
should have T ~ 1. We use as positive direction of the thrust axis the direction of the 
most energetic jet. The thrust axis is defined as the axis along which the projected 
energy flow is maximized. The major is sensitive to the energy flow in the plane 
perpendicular to the event thrust axis, and its direction is defined in the same fashion 
as thrust. The minor axis is the third axis, which is perpendicular to the thrust and 
major axis. The difference between the major and minor values is called the oblateness, 
such that Ob ~ corresponds to an event configuration symmetrical around the thrust 
axis. Figs. ITTk)-c) show the values of the thrust, major and oblateness for the PYTHIA 
6.2 model, after the energy smearing to imitate the LC detector response. The shaded 
band represents the e + e~ — > it — > 6 jets events, scaled according to the predicted cross 
section. It is clear that the values of the thrust and the major have the best sensitivity 
the top-quark production, while the oblateness reflects the properties of the top decays 
in less extent. 

It is mandatory to use an additional information on the invariant mass of the 
decaying system, together with the variables discussed above which mainly focus on 
the topological properties of the event shapes. Since the top mass is approximately 
known, and because the production of two top quarks at s/s > 500 GeV is above the 
top-production energy threshold, an e + e~ — > ti event can be viewed as a simple two- 
body decay with non-overlapping decay products. The top events can be divided in 
two event hemispheres, each of which has an invariant mass close to the top mass. The 
hemispheres can be defined by using the thrust axis, such that particles with positive 
Z component belong the first hemisphere, while particles with Z < form the second 
hemisphere. Fig. fTTH) shows the invariant-mass distribution, M2 = 0.5 (M a + Mb), 
where M a (Mb) is the invariant mass of the hadronic system in the first (second) 
hemisphere. For the e + e~ — > ti — > 6 jets events, the invariant mass has a broad peak 
near the generated (true) top-mass value, while e + e~ events without the all-hadronic 
top decays typically have rather low invariant masses. 
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Furthermore, since the all-hadronic top decay is characterized by a six-jet configu- 
ration, one can use this information for the top selection. Particularly, y™ 1 values of 
the resolution parameter of the Durham algorithm for which an iV-jet configuration 
becomes an N + 1 jet system can also be useful. We have already illustrated that 
the cut on y^ ut alone has a significant impact on the selection of the all-hadronic top 
decays. 

5.1 Standardizing inputs 

The contribution of an input for a NN depends heavily on its variability relative to 
other inputs. Therefore, it is essential to rescale the inputs, so their variability reflects 
their importance. We will rescale some inputs to the interval [0, 1], which is the most 
popular choice for the NN inputs. 

In case of the event shape distributions, we do not use any rescaling, since all these 
variable are already in the interval [0, 1]. For the mass distribution M2, we used the 
function Ai(lnM2 — A2), with A\ = 3.4758 and A2 = 5.0106. This choice was advocated 
in [T^], and here we only adjusted the free parameters A\ and A 2 , such that this new 
variable is zero for m t = 150 GeV, and it equals to unity when m t = 200 GeV. 

The situation with y c ^ 1 values is more complicated, since their values strongly de- 
pend on N. As example, y™* for which 2 jets become a 3-jet system has a significantly 
larger value than y^ 2 - For the all-hadronic top decay, the relative importance of y^f 
values for event configurations with large N is higher than for low N. Therefore, in 
order to treat all inputs on an equal footing, we transformed all y^f values according 
to the expression Y$ = i 2 y^t\ % = 1 . . .6, which increases the relative weights of y^f 
values in the NN input for events with many jets. 

5.2 Training sets 

We have trained the NN to learn just one number in the range between and 1. When 
the NN output is unity, this corresponds to the most probable top candidate in the 
all-hadronic decays. 

A few data sets were used. The first event sample is for the NN training. It consists 
of 40k e + e~ PYTHIA events with the ISR included. One half of this sample contains 
fully inclusive events but without the e + e~ — > ti — > 6 jets process, while another half 
contains the all hadronic decays of ti. The top masses were generated according to 
the Breit-Wigner distribution with the nominal (peak) positions flatly distributed in 
the interval 150 — 200 GeV. The momenta of the final-state particles were smeared as 
discussed in Sect. E2 The second set contains also 40k events, half of them are the fully 
hadronic top decays. This set was used to assess the performance (generalization) of 
the learning, in order to avoid over-fitting. 

5.3 Neural network architecture and the results 

The neural network selection was based on the Stuttgart Neural Network Simulator 
with the logistic activation function 1/(1 + exp(— x)). A feed-forward architecture with 
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one output node representing the probability of the observation of the fully hadronic 
top decays was used. For an input, we use nine nodes representing 1 — T, Mj, the 
rescaling value of the invariant mass, M 2 , and six values of Yj representing the Durham- 
algorithm resolution parameters. We investigated the effect of varying the number of 
hidden layers and nodes in these layers. It was found that one hidden layer with seven 
nodes was sufficient to achieve the best selection performance. 

The neural net was applied to select the top events from the sample discussed in 
Sect. El Events were accepted if the NN output values are above 0.97. Figure IT21 shows 
the NN output, which has a clear peak near unity. 

Figure IT3*1 illustrates the reconstructed top events after the NN event selection. The 
top reconstruction was the same as for FigEl The signal-over-background ratio is larger 
than that without the NN selection, as well as the statistical error on the top mass is 
smaller than for Fig|7||Hl However, a small shift in the reconstructed mass was observed 
which should properly be taken into account after a correction procedure. In addition, 
the shape of the background is less understood than in case without the NN selection. 

6 Accuracy on the top-quark measurements 

The Monte Carlo luminosity used in this analysis corresponds to 16 fb" 1 for y/s = 500 
GeV and 33 fb -1 for y/s = 800 GeV. Our results indicate that the statistical error on 
the top-mass measurement for this luminosity is about 380 — 450 MeV. Thus, for an 
integrated luminosity of 300 fb _1 , corresponding to one to two years of running, the 
proposed method for the all-hadronic decay channel will lead to a statistical error of 
5m t ~ 100 MeV at y/s = 500 GeV. The statistical uncertainty on the mass measure- 
ment for y/s = 800 GeV is by a factor two larger due to a significant missing energy 
from neutrinos. Tables El and El show the results for y/s = 500 GeV and y/s = 800 
GeV. 

The direct measurement of the top-quark width in this decay mode is difficult. The 
reconstructed width is by more than a factor ten larger than the generated (true) top 
width. Therefore, other methods should be investigated to be able to measure the top 
width from the hadronic-final-state signatures of e + e~ collisions. 

One can also roughly estimate the statistical accuracy on the measurements of the 
total cross section for the all-hadronic top-decay channel. This study indicated that 
the number of the top candidates is around 2k for 16 fb -1 . This number of the top 
candidates was obtained from 4224 ti — > 6 jets events, thus the acceptance is ~ 24%. 
For an integrated luminosity of 300 fb -1 , about 80k ti — > 6 jets events are expected 
assuming a t i = 0.6 pb. This leads to ~ 260 ±0.9 fb for the observed cross section, thus 
the relative accuracy is 0.35%. After the event reconstruction, ~ 19200 ± 124 events 
are expected, leading to 64.8 ± 0.4 fb, and to a relative accuracy of 0.6%. 

7 Summary 

For the first time, the all-hadronic top decays were studied using the full TESLA 
detector simulation, after taking into account a realistic contribution from multi-jet 



10 



QCD background. This was done using a new method of the reconstruction, which 
is expected to be less sensitive to experimental and theoretical uncertainties arising 
from the direct measurement of the dijets with the invariant masses close the W mass. 
Note that the method is a very general and can be used for the reconstruction of any 
two particles decaying into six jets in e + e~ collisions. Also, the proposed approach is 
significantly simpler than the top-reconstruction method adopted at the TEVATRON, 
where a kinematic fit is used. The present method takes advantage of the fact that, in 
e + e~ annihilation, the entire hadronic event can be reconstructed. Therefore, essen- 
tially all produced particles can be grouped into jets using an exclusive jet algorithm, 
while energy-momentum conservation provides an additional handle for the top recon- 
struction. In contrast, the reconstruction of top quarks at the TEVATRON and LHC 
is characterized by a large missing energy along the beam directions, and by the use of 
inclusive jet algorithms which combine only a relatively small fraction of the produced 
hadrons into jets. We did not use the double 6-tagging to obtain the top signal; we 
have verified that the fe-tagging reduces the QCD background, but does not have a 
significant impact on the reconstructed width of the top quark. 

The TESLA detector resolution for the top-mass measurement based on the energy- 
flow algorithm is 5 — 5.5 GeV for y/s = 500 GeV. This value is almost independent 
of method used to reconstruct the top quarks. For yfs = 800 GeV, the detector 
resolution is 9.5 GeV. Such an increase is due to higher jet energies, larger overlaps 
between energy deposits used for the energy-flow objects, as well as a larger calorimeter 
leakage. These two numbers can be used for comparisons with other detector designs 
and/or other methods of the event reconstruction. Note that even before the detector 
simulation, a significant width for the three-jet mass spectrum due to fragmentation 
effects is observed, thus the reconstruction of the top width from hadronic jets is very 
unlikely. 

In this paper, we have estimated the typical statistical uncertainties on the top 
mass measurement, as well as the statistical uncertainty anticipated for the top-quark 
production cross section in the all-hadronic decay channel. The method leads to a 
statistical uncertainty of 5m t ~ 100 MeV for the modest value of the luminosity, 300 
fb^ 1 . This uncertainty is compatible with the statistical precision for the most promis- 
ing lepton-plus-jets decay channel at the LHC [2], and is well below the theoretical 
systematic uncertainty anticipated for the pole top mass. An essential aspect of this 
method is to understand limitations arising from experimental systematic uncertain- 
ties. However, this study has not been carried out yet. 
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A E 




A PT 




Am (GeV) 


A P (GeV) 


a) PYTHIA, v/i = 500 GeV 


0.07 


0.04 


0.04 


2 • 10~ 4 


40 


20 


b) PYTHIA, v/i = 800 GeV 


0.4 




0.3 


2 ■ 10" 4 


40 


20 



Table 1: The parameters used for the selection and reconstruction of the top quarks in 
e+e" annihilation at: a) y/s = 500 GeV; h) y/s = 800 GeV. 



PYTHIA, v/i = 500 GeV 


Sm t (GeV) 


BRW width (GeV) 


Gaussian width (GeV) 


Generated 


0.030 


9.5 ±0.4 




a) Reconstructed 


0.105 


9.5 (fixed) 


5.0 ±0.2 


b) Reconstructed (NN) 


0.087 


9.5 (fixed) 


5.2 ±0.1 



Table 2: The statistical uncertainties on the top- mass measurement assuming an in- 
tegrated luminosity of 300 fb" 1 for e + e~ annihilation events at = 500 GeV. The 
table also shows the typical reconstructed widths for the ) without using the 

NN (see Fig.0) and b) with the NN selection (see Fig. ESI). 



PYTHIA, = 800 GeV 


5m t (GeV) 


BRW width (GeV) 


Gaussian width (GeV) 


Generator level 


0.04 


7.7 ±0.4 




Reconstructed level 


0.230 


7.7 (fixed) 


9.5 ±0.4 



Table 3: The statistical uncertainties on the top- mass measurement for e + e~ annihi- 
lation events at s/s = 800 GeV assuming an integrated luminosity of 300 fb" 1 . The 
table also shows the typical reconstructed widths. 
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Figure 1 : The values of ?/g n * for the all-hadronic top decays and for the rest of inclusive 
e + e~ sample (continuum) generated with the PYTHIA model. 
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Figure 2: The distribution of three-jet invariant masses for e + e~ — > tt — »• 6jets events 
at y^s = 500 GeV. The all-hadronic decay channel was generated with the PYTHIA 
model with (open dots) and without (solid dots) the Bose-Einstein effect. The peak 
position and the width were determined from the Breit-Wigner distribution, which gives 
the best y 2 for the fit. A Gaussian convoluted with the Breit-Wigner distribution (with 
the fixed width of 1.39 GeV) is also shown (dashed line). The first-order polynomial 
is used to describe the background. The solid line shows the fit for the invariant mass 
obtained from the PYTHIA models without the BEC32 effect included. 



^300 I i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i i 



> 




250 300 350 400 450 500 550 600 650 

E (GeV) 



Figure 3: The distributions of the total event energies of e + e — > tt events before and 
after the LC detector reconstruction. 
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Figure 4: The reconstructed energies of the final-state leptons for semileptonic top 
decays and for the e + e~ —>■ everything events. 




Figure 5: The invariant-mass distributions for three-jet clusters for e + e~ — > ti process 
at y/s = 500 GeV generated using the PYTHIA model before the LC detector simu- 
lation. The all-hadronic decays were selected and reconstructed using the parameters 
given in Tabled). 
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Figure 6: The invariant-mass distributions for three-jet pairs obtained from the fully 
inclusive e + e~ PYTHIA events. Events e + e~ — > ti were passed through the full detec- 
tor simulation, while the continuum was obtained after the Gaussian smearing of the 
original particle momenta. The top events were selected and reconstructed using the 
parameters given in Tabled), but without the cut on y^ 11 . The fit is based on the 
Breit-Wigner function convoluted with a Gaussian, plus a liner function to describe 
the background. The width of the Breit-Wigner distribution is fixed to 9.5 GeV (see 
Fig. El). 
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Figure 7: Same as Fig. El but after applying the cut y™ 1 > 2 ■ 10' 
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Figure 8: Same as Fig. [7[ but in addition, only three-jet groups are shown which have 
at least one jet pair with | Mjj — M w \< 15 GeV. 
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Figure 9: The invariant-mass distributions for three-jet clusters for e + e~ — > ti process 
at t/s = 800 GeV generated using the PYTHIA model. The all-hadronic decays were 
selected and reconstructed using the parameters given in Tabled). 
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Figure 10: The reconstructed invariant masses for the all-hadronic top-decay channel 
generated with PYTHIA for = 800 GeV after the LC detector simulation. The 
applied cuts as for Fig. [7| The Breit-Wigner width is taken from Fig. cap:tru2. 
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Figure 11: The values of the thrust, major, oblateness and the invariant masses of event 
hemispheres for fully inclusive e + e~ PYTHIA events generated without the tt — > 6jets 
process and for the all-hadronic top decays (shaded aria). 



21 



WW probability 




Figure 12: The output of the neural network for the fully inclusive PYTHIA events 
after the Gaussian smearing of particle momenta to imitate the LC detector response. 




Figure 13: The invariant-mass distribution for three-jet clusters after the NN selection. 
The event selection as for Fig. El 
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